Proximity Effect in Normal Metal - High T c Superconductor Contacts 
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We study the proximity effect in good contacts between normal metals and high T c (d x 2_ y 2- 
wave) superconductors. We present theoretical results for the spatially dependent order parameter 
and local density of states, including effects of impurity scattering in the two sides, s-wave pairing 
interaction in the normal metal side (attractive or repulsive) , as well as subdominant s-wave paring in 
the superconductor side. For the [f 00] orientation, a real combination d + s of the order parameters 
is always found. The spectral signatures of the proximity effect in the normal metal includes a 
suppression of the low-energy density of states and a finite energy peak structure. These features are 
mainly due to the impurity self-energies, which dominate over the effects of induced pair potentials. 
For the [110] orientation, for moderate transparencies, induction of a d + is order parameter on the 
superconductor side, leads to a proximity induced is order parameter also in the normal metal. The 
spectral signatures of this type of proximity effect are potentially useful for probing time-reversal 
symmetry breaking at a [110] interface. 

PACS numbers: 74.45.+c,74.20.Rp 



I. INTRODUCTION 

The proximity effect refers to a broad range of phe- 
nomena related to the leakage of superconducting cor- 
relations into a normal metal in contact with a super- 
conductor. One way of viewing the proximity effect, is 
to consider the decay of the electron pair wave function, 
or off-diagonal retarded Green's function / , from the 
superconductor side to the normal metal side, on the co- 
herence length scalei 

The proximity effect in a normal metal in good con- 
tact with a low-T c superconductor has been studied for 
a long timei*2i2i2i£i£iL& (see also the review in Ref. and 
references therein). On the other hand, the proximity 
effect in a normal metal in good contact with a high-T c 
superconductor has not been considered experimentally 
until very recently. 10-11 

Sharoni et a./J-L performed scanning tunneling spec- 
troscopy on gold coated YBa2Cu 3 07_a (YBCO), and 
found a (pseudo-) gap on the gold side that decayed ex- 
ponentially with distance from the superconductor on a 
scale ~ 30 nm. Kohen et alJ^ measured the conductance 
of a contact between gold and Yi_ x Ca x Ba 2 Cu307_5. 
The experimental results were fitted with BTK-type 12 
(Blondcr-Tinkham-Klapwijk) model generalized to in- 
clude the d x 2_ y 2 symmetry of the order parameter^ as 
well as a subdominant component of the order parameter 
of s or d xy symmetry. The authors of Ref. ^] concluded 
from their fits that the order parameter on the super- 
conductor side might have the time-reversal symmetry 
breaking (TRSB) form d + is. This result was rather 
surprising, since the contacts had predominantly [100] 
orientation and very high transparency, in which case 
TRSB is usually not expected (see e.g. Refs. IHIl5l and 
the reviews in Refs- EHg) . Therefore, the authors spec- 
ulated that maybe there is an unusual proximity effect, 
where the TRSB state is induced by the proximity effect 
at a good contact between a normal metal and d x 2_ y 2 



superconductor with [100] orientation. The authors of 
Ref. ^ later speculated, based on the size of the pseu- 
dogap seen in their STM spectroscopy, that they might 
also have observed this type of new symmetry breaking. 

Motivated by these experiments, we set out to explore 
the proximity effect in normal metal - high-T c supercon- 
ductor structures. We do not find any evidence of a new 
symmetry breaking proximity effect, within our model 
and its parameter space. But, on the other hand, we are 
able to qualitatively explain the pseudogap phenomenon 
seen in the latter experiment in terms of a more conven- 
tional type of proximity effect where the modulation in 
the local density of states is due to the interplay of leak- 
ing superconducting correlations into the normal metal 
and impurity scattering. 



II. MODEL AND METHODS 

We shall describe two types of systems: first, an 
NIS system where a normal metal (N) is coupled to an 
[100] (a-axis) high-T c superconductor (S) facet through 
a high-transparency tunnel barrier (I), as shown in 
Fig-ffl This setup roughly corresponds to the experimen- 
tal situation^ We will also present results for a [110]- 
contact. Second, an INIS system consisting of an a-facet 
with a metal overlayer, shown in Fig. 

Our considerations will be applicable to rather gen- 
eral situations, but in order to match the experimental 
conditions we restrict the parameter space of our model. 
We focus on the low temperature region T = 0.05T C , 
where T c is the superconducting transition temperature 
(T = 4.2K and T c = 90K in the experiment). In the 
S-side we include, besides the d-wave interaction, an at- 
tractive pairing interaction in the s-wave channel with a 
critical temperature T cs — 0.1T C , which is similar in mag- 
nitude to previous estimates from measurementsi2ii2i22i2i 
of TRSB in [110] oriented tunnel junctions of YBCO. We 
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YBa 2 Cu 3 7 . 5 



FIG. 1: We model the experimental setup by a half-infinite 
(a; G [—00, 0]) dirty normal metal in good contact with a [fOO] 
facet of rather clean YBCO (x G [0,oo]). The LDOS in the 
normal metal may be mapped out by STM. 



note, however, that signatures of subdominant pairing is 
not always seen experimentallyi^SiS In the N-side we in- 
clude a pairing interaction that can be either attractive, 
repulsive, or zero. We assume that the gold metal is dirty, 
and include impurity scattering in the self-consistent t- 



matrix approximation 



24.25 



under the assumption that 



the impurities scatter isotropically (s-wave scattering 
only) in the Born limit. This corresponds to the Usadel 
approximation,? 6 which is widely used to describe the 
proximity effect in contacts involving low-T c supercon- 
ductors. Since we are considering the high-T c supercon- 
ductors, which are clean materials with anisotropic pair- 
ing, we can not use Usadel's scheme; instead we solve the 
full Eilenberger equationsSZ^ for the quasiclassical Mat- 
subara matrix Green's function g M (Pf, R; e n ), where p/ 
is the Fermi momentum, R is the spatial coordinate, and 
e n = (2n + \)-kT (n integer) is the Fermion Matsubara 
frequency: 



«e n T 3 - er 



M 



A, 9 
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M 



0, 



(9 M ) 2 



(1) 



Here is T3 the third Pauli matrix in Nambu space, A is 
the superconductor gap (pair potential), and a M is the 
sum of the impurity elastic and spin flip self energies 

~M ,~M 
imp 1 sf' 

Table. [Q 

The gap equation reads 

£„|<W 



& M = ctJ™„ + a^f . We use units H = ku — 1, except in 
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FIG. 2: The second model system is a more conventional 
proximity overlayer, where the LDOS at the normal metal 
surface may be probed by STM. 



where A^ and A s are the interaction strenghs in the d- 
wave and s-wave channels, respectively, in the supercon- 
ductor. The basis function in the d- wave channel is taken 
to have the form r] d (pf) = \/2($%.-f%), J dpfq d (p f ) = 0, 
/dp/|?7ci(p/)| 2 = 1. In the normal metal for attractive 
interaction A„ > 0, while for repulsive interaction A„ < 0. 
For attractive interaction, we can eliminate the cut-off 
u c and the interaction A in favor of the bare bulk super- 
conducting transition temperature T c in the usual way^ 
A" 1 = ln(T/T c )+J] n>1 (n-0.5)- 1 . For repulsive interac- 
tion, this procedure obviously fails, and we have a weak 
cut-off dependence. We take lu c = 30T C , but we have 
checked that our results are not qualitatively changed 
for other values of u> c . 

With the pairing interaction in Eq. (jSJ), the order pa- 
rameter is 



A( P/ ,R) 



A„, N side, 
&dVd(Pf) + A s 



S side, 



(4) 



where all amplitudes are in general complex quanti- 
ties. We choose A^ real, while the phases of A n and 
A s relative to are found self-consistently. We de- 
note the maximum value of the cZ-wave component, 
max{A d r7 d (p / )} = \/2Ad, by A . 

The self energy from elastic impurity scattering has the 

f orm 25 



A( P/1 R)=T £ /^ACp/.p^RJ/^RjO, C(R,e n) = -J-/ d p / ^( P/) R;e„), (5) 



(2) 

where f M is the off-diagonal part of the matrix Green's 
function g M , oj c is an energy cut-off, and J dpf... denotes 
a Fermi surface average. In our model system, we assume 
that the interaction has a simple form 



A„, N side, 



A(p /)P },R) - I \ n d 'vd(Pf)Vd(p' f )+Xs, S side, 



(3) 



and we shall use the mean free path I = VfTi mp as the 
input parameter, which takes the value £ n / s in the N and 
S sides, respectively. We also include spin-flip scattering 
through the self energy^ 

a%(R, e„) = —Jr- / dp f T 3 g M (p f , R; e„)f 3 , (6) 
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where the mean free path l s f = VfT s f is always assumed 
to be large compared to the elastic mean free path. The 
spin-flip process will ultimately set an upper limit on how 
long range the proximity effect isi 

The interface connecting the normal metal side and the 
superconductor side is described by Zaitsev's boundary 
condition^ valid for arbitrary transparency !D (p/) . This 
model describes a wide specularly scattering contact with 
translational invariance along the interface. For simplic- 
ity we consider a two-dimensional system, with quasipar- 
ticles moving within the a6-plane in the superconductor 
which we also denote the xy-plane, with the x-eods ex- 
tending perpendicular to the interface (as in Figs. II 121) 
and the y-axis lies parallel to the interface. We assume 
circular Fermi surfaces in the two sides and take the an- 
gular dependence of the interface barrier as predicted by 
an interface (^-function potential of strength H^, includ- 
ing the possibility of different Fermi velocities in the two 
sides (see e.g. Ref. We assume that the effective 
masses in the two sides are the same (for a discussion 
of effective mass mismatch, see e.g. Ref. fUS). For Zait- 
sev's boundary condition the momentum parallel to the 
interface is concerved, which means that the two different 
trajectory angles n / s of the two sides, measured relative 
to the x-axis, are related as 



Vf s sua ti 8 = Vf n smO n . 



(7) 



We assume that the Fermi velocity is larger in the nor- 
mal metal than in the high-T c superconductor, < a = 
Vf s /vf n < 1. Quasiparticles travelling on trajectories 
with angles n larger than a critical angle 



arcsma, 



(8) 



suffer total reflection. The interface transparency for an- 
gles n < C is 



»(P/) 



Av n V. 



(v n + v 8 ) 2 +4Hf 



(9) 



where v n / s = i>f n / s cosO n / s are the projections on the 
interface normal of the Fermi velocities in the two sides. 
We use the model parameters a and Z = Hb/vf n to pa- 
rameterize the barrier transparency. The transparency 
for perpendicular incidence is denoted So = a/[0.25(l + 
a) 2 + Z 2 }. We concentrate on the high-transparency 
limit, with values given in the experimental paper»ifl 
£>o ~ 0.8-1.0. 

To compute the local density of states, we solve the 
retarded versions of the above equations, which are ob- 
tained by changing the superscript M to R and letting 
ie„ — > e + i0 + in Eqs. (QJ and 0-® ■ Here is i0 + an in- 
finitesimally small positive imaginary number. To speed 
up the convergence of the numerics we have kept this 
imaginary number finite equal to zlO -3 . The gap profiles 
obtained within the Matsubara technique serve as input 
to the calculation of the retarded Green's functions. The 



local density of states (LDOS) is defined as 
A/"(R,e) = J dp/./V(p/,R; e), 
AA( P/ ,R;e) = -^Im {Tr [f 3 g R (p f , R; e)] } 



(10) 



where Aff is the density of states at the Fermi level in 
the normal state. 

The above set of equations for the Green's functions 
and self-energies are solved self-consistently with numer- 
ical methods. This can be done quite efficiently with the 
Riccati parametrization technique ^^S^k 



III. RESULTS FOR THE NIS SYSTEM 

Previously, Ohashi^i made a thorough investigation of 
the proximity effect in a clean system (mean free path 
£ — > oo ) between a normal metal with attractive or repul- 
sive s-wave interaction, in contact with a purely d x 2_ y 2 
superconductor (A s = 0). See also the paper by Bruder4£ 
The central result of Ref. is that for the [100] orien- 
tation an s-wave pair potential is induced in the N-side 
(for A rl 7^ 0), which also leads to spectral changes. For 
the [110] orientation the gap amplitude vanishes. The 
difference between the two orientations is easily under- 
stood by symmetry arguments: the positive and nega- 
tive lobes of the d x 2_ y 2 order parameter contributes to 
the gap equation in the N-side with equal weights (but 
with opposite signs) in the [110] orientation and cancel 
exactly. This cancellation does not occur for the [100] 
orientation, and an s order parameter is induced. How- 
ever, these results does not answer all question raised in 
the experimentsii2*ii Also, the effects of impurities re- 
main unclear. In the following we therefore present an 
extensive investigation of the proximity effect, including 
effects of impurity scattering and subdominant pairing in 
the S-side. 



A. [100] orientation, A n = 

In Fig. |2| we present a representative picture of the 
proximity effect at a [100] contact, in the absence of pair- 
ing interaction in N-side (A„ = 0) but with a subdomi- 
nant interaction channel in the S-side (T cs — O.lTc). In 
part (a) we show the LDOS, and in (b) we show the 
order parameter profiles and some of the characteris- 
tic features in the LDOS. The parameters where chosen 
such that the elastic mean free path in the supercon- 
ductor is £ s — 25£o, corresponding to the typical low- 
temperature residual scattering rate r ~ 1 — 10 ps found 
in experiments^*^ The normal metal side is assumed 
dirty, with £ n = £o = Vf/T c ~ 10 — 100 nm typical for 
metals in general^ and the Au overlayer in the experi- 
ment in particular Here we put vt ~ 10 6 m/s in both 
sides, for simplicity. Note that the superconductor coher- 
ence length £a = Vf/nAo < £o/27r is much smaller than 
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NM 


clean normal metal coherence length £„o = 




1 jim @ IK 




2irfc B T 




normal metal elastic mean free path l n — VfT 


10 — 100 nm 




dirty normal metal coherence length £ n = ^ 


/ hD 


50 — 200 nm © 


1 K 


1 27vk B T 






LTS 


YBCO 


superconductor gap size Ao 




~ 1.5 meV [4] 


20 - 25 meV fl 1.361 


superconductor coherence length £a = 




140 nm 


16 A [37] 


elastic mean free path £ s = W/T 




10 - 100 nm 


250 - 2500 nm [38,39] 


dirty superconductor coherence length £ s — 


1 hD 
V 2A 


27-86 nm 


not relevant (£a <C Is) 



TABLE I: Order of magnitudes of some length scales relevant to the proximity effect in normal metals (NM), low-T c supercon- 
ductors (LTS), and YBCO. The diffusion constant is defined as D = \vjt. For all estimates in NM and LTS we used Vf ~ 10 6 
m/s |40|. The cited LTS gap value is for Nb Q. The maximum d-wave gap amplitude is Ao = V2Ad, and the magnitude given 
here is for the a-axis direction. The short coherence length in YBCO is usually estimated from the upper critical field in terms 
of Ginzburg-Landau theory— H C 2 = <l?o/27r£ 2 , where $o = hc/2e. 
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FIG. 3: The proximity effect at an [100] contact between a 
dirty normal metal (£ n = £o) and a rather clean d x 2_ y 2 su- 
perconductor (l s — 25£o) with a subdominant s-wave com- 
ponent (Tea ~ O.lTc). The pairing interaction in the normal 
metal side is X n — (thus A n = 0). The temperature is 
T = 0.05T C , and the spin flip mean free path is t s f = 100£o 
in both sides. In (a) we show the density of states at x — n£o 
(integer n = —5. ..5) for e G [— 4T C ,4T C ]. There is a sudden 
change at x = because of the interface backscattering for 
So = 0.9 (a = 1, Z = 1/3): the striped filled curve is at 
x — 0~ and the dark filled is at x — + . The density of states 
is normalized to Nf so that it approaches — > 1 for e > T c . 
In (b) we show the spatial dependence of the peak e p (x) in 
the LDOS in the normal metal, the corresponding density of 
states N(tp), the density of states at the Fermi level N(0), 
and the pair potentials A^ and A s in the superconductor. 



the length scale £o = Vf/T c we use. We summarize the 
relevant length scales found in real materials in Table HJ 

If the normal metal was clean, there would not be any 
signatures of the proximity effect in the LDOS^ 1 (since 
the pair potential is zero for A n = 0). In contrast, in the 



dirty case presented here we have clear signatures in the 
form of a pseudogap at low energies, that decays into the 
normal metal on the thermal coherence lenght scale = 
\J D I (2ttT), where D = Vf£ n /3 is the diffusion constant 
in the normal metal. To quantify the effect we show in 
Fig. ISIb) the spatial dependences of the pseudogap peak 
position e p (x) (solid black line), the LDOS at the peak 
position Af(e p ) (red dashed-dotted line), and the LDOS 
at zero energy Af(0) (green dashed-double-dotted line), 
where it is suppressed the most. The maximum deviation 
from the normal metal density of states is about 25%, 
which is much smaller than the 100% effect predicted in 
low-T c s-wave superconductors^iiSiLSiS at zero energy. 
This difference is due to the nodes in the d x i_ y i -wave 
order parameter and the cancellation effects between the 
positive and negative lobes, which substantially reduce 
the influence of the unconventional superconductor on 
the LDOS in the normal metal. 

The presence of a proximity effect at the junction sig- 
nals a non-zero Fermi-surface average 

|dp / / fl (p / ,R;e)^0 (11) 

in this region. This average is non-zero near the inter- 
face in both sides. Its non-zero value is intimately re- 
lated to the suppression of the d x 2_ y 2 component , 41 i 43 i 44 
Since / fl (p/,R;e) cx %(P/); a non-zero Fermi surface 
average can only be obtained for a spatially dependent 
d 2 .2_ !/ 2-wave gap in a self-consistent solution where non- 
locality is taken into account. The average is also en- 
hanced when the back-scattering probability at the inter- 
face 9^o = 1 — £>o is non-zero, since then the importance 
of high angle trajectories > 7r/4 (for which the <i-wave 
order parameter is negative) is reduced compared to the 
low angle trajectories < 7r/4 (for which the order param- 
eter is positive). 

In Fig. we couple back the non-zero Fermi surface 
average of the off-diagonal Green's function in the su- 
perconductor side by having a small subdominant pair- 
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FIG. 4: The same as in Fig. [3] but with repulsive pairing 
interaction in the normal metal side A n = —5, which result in 
an induced pair-potential A n in the normal metal. 

ing interaction (X s > 0). Wc then obtain a non-zero 
subdominant order parameter component A s near the 
interface in the superconductor. However, this compo- 
nent is always in phase with the dominant <i-wave, and 
time-reversal symmetry is always conserved at the [100]- 
contact in the parameter space we have explored (mean 
free paths ranging from the clean limit £ n — > oo to the 
dirty limit l n <C £ n o with a wide variety of barrier trans- 
parencies and subdominant pairing interactions). The 
effects of the s-wave gap on the local density of states is 
marginal. 

B. [100] orientation, A n < 

In Fig. 0] we present the proximity effect at a [100] 
contact, for the case of a repulsive pairing interaction in 
the N-side. We have in mind that the counter electrode in 
the experiment ™ was Au which is not a superconductor 
and might have a repulsive effective pairing interaction. 
The d x 2_y2-wewe pairing correlations leaks over to the N- 
side, which leads to the non-zero Fermi-surface average 
in Eq. i|ll|) . That average is coupled back though the 
repulsive pairing interaction A„, which leads to a non- 
zero pair potential A„, that decays into the bulk normal 
metal. The subdominant component A s on the S-side is 
not qualitatively affected by the induced pair potential 
on the N-side. The opposite is also true, the induced 
A„ is not qualitatively affected by the presence (A s ^ 0) 
or absence (A s = 0) of the component A s . Thus, the 
proximity induced pair potential on the normal metal 
side is mainly due to the dominant d x i_ y i component. 
We also note that the real combination d + s is always 
favored in the system, as in the case A„ = above. 

As for the LDOS, the induction of the gap amplitude 
A„ changes the density of states near the contact. There 
are new low-energy states, which are formed because of 
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FIG. 5: Proximity effect at a [100] contact with a normal 
metal with attractive pairing interaction (T n = 0.0 1T C < T = 
0.05T C ), for high barrier transparencies Do = 1 (a = 1, Z = 0) 
in (a) and Do = 0.9 (a = 1, Z = 1/3) in (b). All other model 
parameters are the same as in Figs. [3] Only the density of 
states in the normal metal at x £ [— £o,0] (in steps of 0.1£o) 
for e G [-4T C ,4T C ] is shown in (a)-(b). The solid (dashed) 
lines in (c) refer to the contact with So = 1 (So = 0.9). 

the sign change of the order parameter field for quasipar- 
ticles travelling on low-angle trajectories < tt/4 which 
connect the negative gap A„ on the N-side with the pos- 
itive lobe of the d x 2_ y 2 order parameter on the S-side. 
This is analogous to the zero-energy bound states formed 
at a [110] surface of a d x 2_ y 2 superconducor^ and was 
also discussed in the study of the superclean case^i In 
our case, the impurity scattering broadens the states con- 
siderably. Thus, compared with the case A„ = above, 
the LDOS peak position e p (x) is drawn closer to the 
Fermi surface. At the same time, the overall modulation 
of the LDOS in the normal metal is reduced compared to 
the case without the pairing interaction. Thus, contrary 
to our first expectations, an induced non-zero pair poten- 
tial A„ in the normal metal reduces the proximity effect 
in the LDOS, as compared to the case A„ = 0. This is a 
unique feature of the dirty system. 



C. [100] orientation, \ n > 

For the case of attractive interaction in the N-side, 
A n is positive which results in a sign change of the or- 
der parameter A n compared to the repulsive interaction 
case considered above (see Fig. EJ. At the same time 
the order parameter amplitude becomes larger and de- 
cays slower compared to the repulsive interaction case 
(analogous results were also found for the clean case in 
Ref . l41|) . As for the LDOS, the weight of the low-energy 
states is larger because of the larger pair potential, al- 
though the trajectories contributing to these states are 
at high angles > 7r/4. For a fully transparent interface, 
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the states are centered at zero energy, Fig.[5fa). For non- 
zero back-scattering at the junction, the states are split, 
see Fig.EJb), where the split size depends on the interface 
transparency (£>o = 0.9 in the figure). The total effect of 
the proximity effect is enhanced for the case of attractive 
paring interaction in the normal metal, compared to the 
case without interaction in Fig. [3J The LDOS variations 
can be as high as 80% for our choice of parameters, see 
Fig. He). 

Given the fact that there are low-energy states we are 
tempted to think that a TRSB combination d x 2_ y 2 + is 
could be formed (at least for Do = 1 f° r which the split- 
effect is absent), instead of the real d x 2_ y 2 + s combi- 
nation found above. If this was the case, the energy of 
the system might be lowered, similar to what happens at 
[110] surfaces or Josephson j unctions However, 
for the proximity contact we were never able to stabi- 
lize a TRSB state within the parameter space we have 
considered. We understand this result as due to a phase 
locking between the s-wave gaps in the N and S sides, to 
the dominant d x 2_ y 2, that prevents a relative phase 7r/2 
between any two components to be formed. Compared to 
the Josephson junctions case, the two sides of the prox- 
imity contact are not independent: there is no bulk order 
parameter on the N-side that makes it possible to estab- 
lish a phase difference 7r/2 over the junction. As for the 
subdominant s-wave order parameter on the S-side, it is 
due to the spatially varying d x 2_ y 2 order parameter, and 
is not induced by the bound state as in the [110] surface 
case. We also note that spontaneous symmetry break- 
ing in the form of a paramagnetic instability 34 - 48 - 50 51 - 52 
is suppressed by impurity scattering and normal back- 
scattering at the junction. 



D. [110] orientation 

We now turn to the [110] orientation. In previous stud- 
ies of the clean system^ no gap amplitude was found in 
the normal metal side for reasons of symmetry. As for 
the LDOS in the N side, it is always unchanged in a clean 
system when the gap amplitude is zero. For A„ = 0, we 
find no signatures of the proximity effect in the LDOS in 
the dirty case either, which also follows from symmetry 
arguments. However, when we include a subdominant 
component in the superconductor side at a contact with 
[110] orientation, it will be formed tt/2 out of phase with 
the dominant, and we have a TRSB stat o 14 i 15 (this is 
opposite to the [100] case above). We note that some 
backscattering is necessary at the junction, i.e. if the 
transparency Do = 1 there are no zero-energy states and 
no TRSB. As for the proximity effect in the N-side, high 
transparency is favorable. Thus, we concentrate on a 
window of intermediate transparencies, in order to es- 
tablish a proximity effect between a TRSB d x 2_ y i + is 
state at a [110] interface to a normal metal, see Fig. [5] 
Clearly, the established is-part in the S-side, leads to a 
proximity induced is component also in the N-side that 
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FIG. 6: The proximity effect at a [110] contact with T cs = 
0.1T C . A subdominant component developes near the surface 
in the TRSB combination d + is at the S-side. For inter- 
mediate transparencies, here Do = 0.3 (a = 1, Z 2 = 7/3), 
the subdominant component induces an is-wave pair poten- 
tial at the N-side (A„ = -5). We have T = 0.05T C , £ n = £ , 
i s = 25£o, and t s} = 100&. 



decays into the bulk normal metal, although the effect 
is rather small and short range. The details of the win- 
dow of transparencies depend on the mean free paths, as 
well as the interaction strengths. Typical values where 
there are TRSB and at the same time visible effects in 
the LDOS on the N-side are Do ~ 0.1 — 0.5, for our choice 
of parameters. The double peak structure on the S-side, 
see the (a)-part of the figure, decays on the coherence 
length scale, together with the subdominant component, 
red-dashed line in (c). Both this peak structure, and the 
peak structure in the LDOS at the N-side, [black dashed 
line in (c)], are related to, but not direct images of, the 
s-wave pair potentials «A S and iA n because of the im- 
portance of the impurity self-energy. 

These results have the potential to serve as a probe 
of TRSB at [110] interfaces: as the temperature is de- 
creased, the spectrum in the normal metal overlayer re- 
mains unchanged (equal to Mf for all energies) until the 
TRSB transition is reached, below which a small pseudo- 
gap develops as shown in Fig. |SJ This probe is however 
not very quantitative and will not give direct information 
about the size of A s , because of the masking effect of the 
impurity scattering. Rather, it would serve as a yes/no 
experiment, as to the existence or absence of TRSB at 
low temperatures. 

The drawback of the above considerations is that they 
rely heavily on symmetry arguments. Thus, to a much 
higher degree than other results in this paper, non- 
specular interface scattering 53 - 54 - 55 - 56 - 57 - 58 could be detri- 
mental. 

Another caveat is the possibility of d-wave pairing cor- 
relations (attractive or repulsive) in the normal metal. 
An induced d-wave pair potential leads to spectral 
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FIG. 7: Quasiparticles travelling on high-angle trajectories 
On > dc suffer total reflection at the NIS interface because 
of the Fermi velocity mismatch Vf n > Vf s (the trajectory to 
the right). On other trajectories 6 n < C , quasiparticles are 
partially reflected in accordance with Eq. ©. 



changes in the normal metal for the [110] orientation also 
in the absence of TRSB, thereby spoiling the symmetry 
argument. The LDOS is enhanced at low energy by the 
formation of low-energy Andreev states broadened by im- 
purity scattering, if the ri-wave pairing correlations are 
oriented as [110] relative to the interface. This is similar 
to what happens at Josephson junctions between two d- 
wave superconductors^ However, it seems unlikely that 
spatially homogeneous ri-wave pairing correlations with 
a well-defined orientation relative to the interface could 
exist in a normal metal such as Au. 

However, these complications are temperature- 
independent and should not appear suddenly at the low- 
temperature TRSB transition. Thus, a smoth but rather 
sudden change in the normal metal spectrum near the 
interface as a function of temperature can be a signal of 
TRSB. 



IV. RESULTS FOR THE INIS SYSTEM 

Let us now study the change in the density of states 
in a normal metal overlayer in good contact with a [100] 
facet, see Fig. EI Because of the absence of a bulk bound- 
ary condition in the normal metal overlayer, we solve 
Eq. iteratively until self-consistency is reached in the 
overlayer, under the asumption of specular reflection at 
the vacuum- normal metal surface (at x = —d). Includ- 
ing the possibility of a Fermi velocity mismatch, we have 
two types of trajectories when we solve Eq. JIJ depend- 
ing on if the trajectory angle 9 n is larger or smaller than 
the critical angle for total inner reflection in Eq. (JSJ) . We 
show the trajectories in Fig. [3 Note that when we draw 
the trajecories in Fig. we do not imply that quasiparti- 
cles move ballistically. Rather, the lines denote trajecto- 
ries along which we solve Eq. QJ. Quasiparticles diffuse, 
which is described by solving the impurity self-energies 
in Eqs.JSj-jni sclf-consistently with the Green's function 
in Eq. (JTJ. 

Because of the severe isotropization effect caused by 
impurity scattering, we do not expect non-specular scat- 
tering to drastically alter our results. In fact, similar 
results for the LDOS in the superconductor as we find 
was already reported by Golubov and Kupriyanov^&M 




4./T 



FIG. 8: The proximity effect at a [100] facet with a d = 2£n 
thick normal metal overlayer. We have T — 0.05T C , T cs = 
0.1T C , A„ = -5, i n = Vfn/Tc = 2v fs /T c = 2£o, t 3 = 25£o. 
The Fermi velocity mismatch is «/ n = 2vf s (Z — 0), which 
leads to £>(#„ = 0) w 0.9 and C « 30°. 



They used the dirty normal metal overlayer as a model 
for a diffusive boundary condition for the superconduc- 
tor. In their case, however, the dirty layer thickness was 
assumed small, d -C V£o^n- Here we relax this assump- 
tion and study in more details the effect of the supercon- 
ductor in the normal metal. 

We present typical results for the order parameters and 
the LDOS for a [100] contact in Fig. We find that the 
LDOS is almost independent of coordinates in the nor- 
mal metal, but depends on the thickness of the overlayer. 
There are quasibound states, seen as peaks in the LDOS, 
reflecting interference between electrons and holes in the 
dirty normal metal in contact with the superconductor. 17 
The energies of these states depends on the thickness of 
the overlayer. In a normal metal in good contact with 
an s-wave superconductor, there is a minigap E g of the 
order of the Thouless energy E t h = D/d 2 (D = \vf£), 
below which the density of states vanish (for a discussion 
see e.g. Ref. and references therein). Because of the 
nodes of the d x 2_ y 2 order parameter and impurity scat- 
tering, this effect is removed and only a suppression of 
the low-energy density of states remains. As in the NIS 
case discussed in the previous section, the LDOS in the 
normal metal mainly reflects the influence of the impu- 
rity self-energy and does not contain direct information 
about the induced pair potential A„ for non-vanishing 
interaction A„ ^ 0. 

In Fig. El we show typical results for the [110] contact. 
The peaks in the LDOS in the overlayer are shifted?! 7 , 
as compared with the [100] contact. This reflects the 
formation of the zero-energy statesii 7 ^ For lower trans- 
parencies, the TRSB state can be favorable, which leads 
to a split of the zero-energy peak, similar to the situation 
in Fig. EJ 

In Fig. ^3 we plot the surface LDOS and the conduc- 
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resolved surface local density of states Af(x = — d, e) 




FIG. 9: The same as in Fig- |H1 but for the [110] orientation. 
For these model parameters only a d-wave order parameter 
exist. The peak positions in the normal metal are at e p = 
{0,0.96T C ,1.98T C }. 




FIG. 10: (a) The surface local density of states and (b) the 
conductance as it would be measured by STM according to 
Eq. H12H for a #* lp = 10° wide tunnel cone, for several dif- 
ferent thicknesses of the normal metal overlayer ranging from 
d — £o to d — 4£o- The mean free path is £ n — 2£o and total 
isotropization by the impurities has not yet occured which 
leads to an enhancement of the modulation in the surface 
density of states as seen in the STM conductance [note the 
different scales in (a) and (b)]. The arrows indicate how the 
peak positions move toward the Fermi level for increased over- 
layer thickness. The other model parameters are the same as 
in Fig. |H| Note that the red dashed curves corresponds to the 
data in Fig. [H] 



tance at the [100] contact as it would be measured by an 
STM tip, as a function of layer thickness. For simplicity 
we model the tunnel barrier between the tip and the over- 
layer by a tunnel cone of width 29 l ^ nel with transparency 
J)q P for angles within the cone and transparency for 
angles outside. The conductance is related to the angle 



-00 = — 

dV y ' f$» 



d9 n D 



tip 



N{-d, 9 n ;eV) 



a/; 



(12) 



where i?^ lp is the resistance of the tip-sample contact in 
the normal state (we normalize J d9 n S) tlp (9 n ) = 1). For 
intermediate overlayer thicknesses, total isotropisation of 
the LDOS by the impurities might not have occured, and 
a small tunnel cone can enhance the structures in the 
LDOS as seen by the STM, since trajectories close to the 
interface normal is selected. 



V. SUMMARY AND CONCLUSIONS 

We have studied the proximity effect in good contacts 
between d x 2_ y 2 superconductors and normal metals. We 
have investigated the effects of impurity scattering, a sub- 
dominant component of s-wave symmetry in the super- 
conductor side, and a proximity induced order parameter 
in the N side. For the [100] orientation, we always find 
a real combination d x 2_ y 2 + ,s in the superconductor. At 
the same time, the N side gap amplitude is phase locked 
to the d-wave order parameter on the S-side. On the 
other hand, for a [110] contact with intermediate values 
of the transparency, a TRSB d x 2_ y 2 +is order parameter 
on the S side, induces an is gap amplitude also on the 
N side. We note that the subdominant is component in 
the S side is not due to the proximity effect; rather, it 
is due to the formation of the bound state at the [110] 
interface in the usual wayi^i^ The signatures of the prox- 
imity effect at [110] contacts could serve as fingerprints 
of TRSB, with the caveat of possible breaking of the "se- 
lection rule" by e.g. non-specular scattering or d-wave 
pairing correlations in the normal metal. 

We finish with a comparison with the experiments^^ 
We can not find evidence of the unusual symmetry break- 
ing proximity effect at [100] contacts that was proposed 
m Ref.Hol The proximity induced pair potential in the 
normal metal is due to the suppression of the d-wave or- 
der parameter and no non-trivial relative phase develops. 
However, a serious comparison with the experiment of 
Kohen et alML can not be done since that was a transport 
experiment where the conductance was measured. At the 
same time, the analysis 1 — in terms of a BTK-type theory 
with postulated pair potentials and no considerations of 
real non-equilibrium effects, rather than self-consistently 
computed quantities, might be too simplistic. Especially 
when the contact has high transparency. For example, 
a dirty normal metal (instead of a clean metal as in the 
BTK theory) in contact with a d-wave superconductor 
can lead to quite different results for the conductance, 
as shown theoretically recentlyi^SS Those considerations 
were limited to point contacts with the superconductor 
a purely d-wave (A s = 0), super-clean (i s — > oo), unper- 
turbed reservoir. 
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On the other hand, our results is in qualitative agree- 
ment with the experimental results of Sharoni et al. 11 
The modulations of the density of states consists of a 
suppression of the density of states at the Fermi level 
(but not to zero), with a quite well defined peak posi- 
tion e p . These structures decay into the normal metal on 
the scale £ n . However, an unambiguous quantitative fit 
between theory and data is hard to make for several rea- 
sons. On the theory side, even within our simple model, 
the local density of states is sensitive to many param- 
eters, such as the mean free paths, the interface resis- 
tance, the Fermi velocities, as well as pairing interaction 
strengths. The choice of parameters is not unique. On 
the experimental side, only a limited number of spectra 
were shown, although many points e p (x) were reported. 
In particular, the spectra for small distances were not 
shown, for which e p were reported to be spectacularly 
large ~ 15 meV ~ 0.75Ao, which is hard to reproduce 
in theory (e p ~ 0.25A in Figs. I3I6|) . We note, how- 
ever, that similar (but not as large) discrepancies between 
theory and experiments were recently reported^ also for 
low-T c contacts (Nb-Au). Although in that case, the dis- 
crepancy was the largest for large distances. Another 
problem is that the modulation in the LDOS appears to 
be larger than we find (unless the effective pairing in- 
teraction in Au is attractive which seems unlikely). In 
particular J\f(e = 0) is suppressed more than we find for 



the NIS system. On the other hand, the experimental 
situation was not clean cut: a mixture of the NIS and 
INIS systems appears to be relevant. For large (small) 
distances from the a-facet, the setup appears to be NIS- 
like (INIS-like) . It is also unclear where along the ir-axis 
in our Fig. ^ the spectra were actually recorded. From 
their Figs. 1-2 it appears that the spectra were actually 
taken at points corresponding to x > in our Fig. ^ i.e. 
into the YBCO c-axis. 

Despite these complications, it is quite plausible that 
the experimental results could be explained by the more 
conventional type of proximity effect we have discussed 
in this paper, rather than an exotic proximity effect that 
involves symmetry breaking. 
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